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ABSTRACT 

We analyze the temporal evolution of accretion onto rotating black holes subject to large-scale magnetic 
torques. Wind torques alone drive a disk towards collapse in a finite time ~ tffEk/EB, where tff is 
the initial free-fall time and / Eb is the ratio of kinetic-to-poloidal-magnetic energy. Additional spin- 
up torques from a rapidly rotating black hole can arrest the disk's inflow. We associate short/long 
gamma-ray bursts with hyperaccretion/suspended-accretion onto slowly/rapidly spinning black holes. 
This model predicts afterglow emission from short bursts, and may be tested by HETE-II. 

Subject headings: black hole physics — gamma-rays: bursts and theory 



1. INTRODUCTION 

Active black holes are believed to be surrounded by mag- 
netized disks or tori. In particular, such systems are the 
most favorable inner engines of cosmological gamma-ray 
bursts. The BATSE catalogue shows a broad, bimodal dis- 
tribution of GRB durations (Tgo) centered on short bursts 
of ^ 0.3 s and long bursts of '--^ 30 s (Kouveliotou et al. 
1993; Paciesas et al. 1999). The cosmological origin of 
GRBs and their rapid variability indicate an association 
with stellar-mass progenitors (Woosley 1993; Paczyhski 
1991, 1997, 1998; Fryer et al. 1999). Proposed scenarios 
include coalescence of compact object binaries (Paczyhski 
1991), and core-collapse of massive stars (Woosley 1993). 
Both scenarios involve a black hole-disk or -torus sys- 
tem. This disk will be magnetized, due to seeding from 
the progenitor's magnetic flux and disk dynamo amplifica- 
tion (Narayan, Paczyhski, & Piran 1992; Paczyhski 1997; 
Brandenburg et al. 1995; Hawley, Gammie, & Balbus 1996; 
Hawley 2000). Various lines of observational evidence sup- 
port the coUapsar scenario, at least for long bursts (e.g. 
van Putten (2001)); here, we shall suggest that both short 
and long bursts can be explained in a unified model with 
durations extended by coupling to a high-spin black hole. 

Accretion provides a definite power source for GRBs, 
but the black hole's spin energy, e.g., 2Mqc^ for a 7Mq 
black hole, can be at least as important. Efficient spin 
energy extraction requires an external magnetic field sup- 
ported by surrounding matter. Magnetically mediated 
energy extraction mechanisms include radiation in DC 
Alfven waves (Blandford & Znajek (1977); Thorne et al. 
(1986); see also e.g. Lee et al. (1999); Beskin & Kuznetsova 
(2000); Krolik (2000)), or superradiant scattering of fast 
magnetosonic waves (Uchida 1997; van Putten 1999). 

The evolution of the disk or torus surrounding a black 
hole depends crucially on magnetic torques. Negative 
torques are expected from magnetic winds, which can 
dramatically stimulate accretion rates (e.g., Blandford 
& Payne (1982); Lee & Kim (2000)). Positive torques 
T = —Jh may derive from coupling to the black hole's an- 



gular momentum Jh by equivalence in poloidal topology 
to pulsar magnetospheres (van Putten (1999); Brown et 
al. (2001)), and may be key to driving hypernovae (Brown 
et al. 2000). The disk may also redistribute its mass and 
angular momentum via internal, turbulent, magnetic and 
Reynolds stresses (Balbus & Hawley 1998). The relative 
contributions from external (wind/black hole) stresses vs. 
internal stresses is presently unknown, and may depend 
crucially on the net poloidal flux introduced to the disk 
(see e.g. Matsumoto (1999)). 

In this work, we focus on the duration of magnetized ac- 
cretion regulated by external torques. We show that wind 
torques alone drive the disk to a finite-time collapse singu- 
larity. The additional positive torque from the black hole, 
when its angular velocity exceeds a critical value, may ar- 
rest accretion onto the black hole for the duration of its 
spin-down lifetime. This introduces two cases: a short- 
lived state of "hyperaccretion" for low-spin black holes, 
and a long-lived state of "suspended-accretion" for high- 
spin black holes. The respective time-scales are consistent 
with the bimodal distribution in GRB durations for plausi- 
ble values of two presently-uncertain ratios: the kinetic-to- 
poloidal-magnetic energy in the disk and the black hole-to- 
disk mass. Numerical simulations soon within reach may 
establish the likely ranges of these ratios. 

Previously, bimodality in GRB durations has been at- 
tributed to, e.g., hydrodynamic time-scales in Newtonian 
vs. relativistic reverse shocks in shells decelerating in the 
ISM (Sari & Piran 1995), or formation of millisecond pul- 
sars spinning above or below a gravitationally unstable 
limit (Yi & Blackman 1998). Previous studies of Poynt- 
ing dominated winds in GRBs have focused on plasma 
processes and possible connections to gamma-ray emission 
(e.g. Usov (1992, 1994); Meszaros & Rees (1997);Black- 
man & Yi (1998);Lyutikov & Blackman (2001)). 

2. ORBITAL EVOLUTION OF ROTATING, MAGNETIZED 
RINGS SUBJECT TO WIND TORQUES 

Magnetic disk-winds transport angular momentum by 
Reynolds and Maxwell stresses, which for outgoing (in- 
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going) flows have the same (opposite) sign. A force-free 
limit is defined by neglecting Reynolds stresses. Applied 
to outgoing winds, Maxwell stresses —RBrB^/An provide 
a lower limit on the disk's accretion rate, where B de- 
notes the magnetic field with components in spherical co- 
ordinates {r,9,(j)), and R = rsinO. The force- free limit 
represents pure "magnetic braking." 

A magnetized disk ring, with mass SM and radius TZ 
in quasi-Keplerian rotation at rate fl = {GMh /'R?Y^'^ , 
evolves as it sheds angular momentum in a wind. Assum- 
ing radially asymptotic flow, the magnetic torque 5tb on 
the ring equals (Goldreich & Julian 1969): 



5tb = - sin^ e [d^^A^) 5A4 



(1) 



Here, fi — cos 0, the poloidal magnetic field is Br = 
—r~'^d^A^, and 27r(5A^ is the magnetic flux through the 
ring. The toroidal field is B^ = —flRBr/c (see the 
wind asymptotics of Goldreich & Julian (1969); Michel 
(1969) for aligned rotators; for disk winds, see Blandford 
(1976); Lee & Kim (2000)). The ring evolution equation, 
6tb = {l/2)fimZ6M, yields 



nn = - sin^ e (a^A^) 



(SAis 
\5M )' 



(2) 



The evolution of a collapsing ring depends on its de- 
gree of magnetization and the magnetic wind's geome- 
try. A split-monopolc geometry (SMG) for the asymptotic 
poloidal field is a good approximation when [Srl/I-B^l >> 
1, true at the onset of accretion for modest rotation rates. 
As the disk collapses and spins up, however, magnetic field 
line winding produces large toroidal magnetic fields, and 
a force-frce-toroidal-field geometry (TFG) becomes more 
appropriate. SMG/TFG represent two limiting cases, for 
which simple solutions to equation (2) may be obtained. ^ 

Solution for split-monopole geometry — In SMG, we 
have A^ = A{1 — /u//xo), with 9o = arccos(y^o) the angle of 
the innermost streamline with respect to the pole; the out- 
ermost streamline lies on the equator. With w = TZ/TZq, 
for TZo the the initial ring radius, the solution of equation 
(2) describes a finite-time singularity: 

zu{t;A^)=(l-t/tJ^'y^\ (3) 

where = {1/4){5J° /\St%\), and 5J° = SMfloTll and 
5t^ denote the initial ring angular momentum and torque. 

The initial ring kinetic energy is 6E^ = 0,^5 / 2. The 
energy in the magnetic field of a (double) wedge-shaped 
portion 5A^ rooted at TZa is 5E^ = (1/2) / B^d^r^dr = 
M^|a^A^|/(27^o) = |^TB|c(27^o^^o)"'(sin^?)-2. For split- 
monopole geometry, 6E% 6E% ^ = A6A4/{2no). The 

initial freefall time is t// = 7r/(rio2\/2)- For Rg = 
2GMh/c^ the Schwarzschild radius of the black hole, the 
collapse time in SMG becomes 



,smg 



SEk _ tff/TT 



2cRs sin^ e 6Eb,s 
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7^o^'^' 5Ek 
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B,s 



(4) 



In the above, superscripts on the energies are dropped 
since SE°/SE%^^ = {5M / 5A^{GMh I A) is an evolution- 
ary invariant. Equation (4) shows that accretion is more 

^ Although, as shown by Blandford & Payne (1982) for (baryonic) 
poloidal fields, cylindrically-coUimated MHD winds have been found 



efficient with increasing sm.0; outer-disk rings have mag- 
netic field lines closer to the equator than inner-disk rings, 
providing greater moment arms for the magnetic torque. 

Solution for toroidal field geometry — Under large mag- 
netic hoop stresses, magnetic flux surfaces adjust their lat- 
itudes such that the toroidal magnetic field in the wind 
approaches a force-free configuration with B^ oc R~^ (see 
e.g. the protostellar disk wind in Shu et al. (1995)). The 
magnetosphere is then effectively current- free except along 
the pole and the equator, with the poloidal pole current 
Ip producing —2Ip/{cR) elsewhere. Hence, d/^A^ = 

-r'^Br = -2Ip[nsm^ 9]-'^ so that \Stb\ = 2IpSA^/c. 
Here, Ip relates to the total normalized fiux A = J dfi df, A0 
by Ip ~ —A^l(TZin)/{2\ ln(f?o/2)|), where TZm is the inner- 
most disk radius. 

From equation (2), the TFG limit gives 07^7^ = 
-{4:Ip/c){SA^/6M). Using ^^(7^) = ^7(7^o)ro-3/2 we find 
the finite-time singularity solution: 

w{t;A^) = {l-t/ty''f, (5) 

where 

<5J° dEk 



^tfg 



B,s 



Here, Jq = {c /V2)7lo^ J dOiQ sin is ~ 

(7^f„/i^57^g)^^^ I ln(6'o/2)|. Going from 7^,„ to the outside 
of the disk at TZd, the radial factor in Xq decreases from 
(7^,:„/7?,,)l/2 (which is order-unity) to {nlJRsTllf/'^ 
(which is of order TZin/TZd)', the product tffl^ varies 



oc 7^J/^ Note that 5E% in TFG is larger than 5E% ,^ by a 
factor \d^A^\IA ~ ^^(7^i„)[^](7^o) sin2 e\ ln(0o/2)|]-^. 

For the inner disk, t^J^ « tj"^ due to the geomet- 
rical factor 46'qI ln(6'o/2)| << 1; the enhanced efficiency 
arises from squeezing of the poloidal field lines toward 
the pole by toroidal hoop stresses. For the outer disk, 
^tfg . ^sn^g ^ 4(7^„,/7^rf)3/2 1 ln(0o/2) | , which is of Or- 
der unity or less unless the effective beaming angle 6*0 is 
extremely small. Thus, although inner-disk accretion is 
greatly enhanced when the magnetic wind becomes more 
stratified toward the poles, the overall lifetime set by 
outer disk accretion is not as sensitive to these geomet- 
ric changes. 

Dimensional estimates of accretion times — Dimension- 
ally, the minimum inner-disk accretion time, from TFG, 
is 

' (_Mh SEk 



tin = 0.011 sgiOo) 



(7) 



^100 kmy \10MqJ 6Eb,, 

where .g(e'o) = (1 + ln(6'o/l°)/ ln(0.009)) and Oq is the an- 
gle of the innermost wind field line. Correspondingly, the 
maximum outer-disk accretion time is 

*. = 0.057s(,^)7^V'#^, (8) 



1000 km 



SE 



B,s 



assuming TZd » TZin with SMG; otherwise, if TZd ~ TZin, 
we would use the TFG result, which increases the time by 
a factor ^ 4| ln(^o/2)| consistent with equation (7). These 
limits will bracket the interval over which disk accretion 
takes place. 

MHD disk winds, field line focusing may produce nonradial asymptotic 
to have relatively low asymptotic flow speeds (Ostriker 1997). 
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3. DISK ACCRETION AND BLACK HOLE SPINDOWN TIMES 
COMPARED 

A ring spends the greatest portion of its lifetime near 
its initial radius; during the subsequent rapid infall phase, 
magnetic flux is approximately conserved. For fiducial 
estimates, we shall normalize to an initial kinetic-to- 
poloidal magnetic energy ratio ~ 100 for all rings, as found 
from Hawley (2000) 's numerical simulations of global, hot, 
MHD accretion disks. This energy ratio certainly varies in 
real systems; estimates based on simulations with differ- 
ent physics or initial conditions (e.g. heating from resis- 
tive dissipation, as in Stone & Pringle (2000); net poloidal 
flux, as in Hawley, Gammie, & Balbus (1996)) may de- 
crease this ratio toward unity. The natural diversity in 
mean magnetizations may in part be responsible for the 
broad distributions of the two duration classes of GRBs. 

Time-scale for short GRBs — Applying 5Ek/ 5Eb,s = 
10 — 100 in equation (7) or (8) (for a disk of size 
;S 1000 km), the accretion time will be ^ a few seconds - 
comparable to that of the 2s event GRB 000310C (Jensen 
et al. 2001). In this scenario, short bursts represent hy- 
peraccretion onto a slowly-spinning black hole. 

Time-scale for long GRBs — What accounts for the 
long-duration bursts? If the black hole is initially 
rapidly spinning, then the burst may persist for the time 
tH required for angular momentum extraction by the 
Blandford-Znajek (1977) and related processes. Given a 
horizon flux the black hole torque for open flelds 

is Ch ^ Q.hA\/c (sec Thornc et al. (1986)). To esti- 
mate AH^ assume that the black hole field strength is 
comparable to that from the inner part of the disk. In 
TFG, we find Ah « A(2| ln(6»o/2)|)-i so that Ch ^ 
0//A2(4c| ln(6'o/2)p)"^ For comparison, the total disk- 
wind angular momentum luminosity is equal to r_ = 
/drs; i.e. C'J'> = a{n,^)A^{c\\n{6a/i)\y\ The spin- 
down time tn = 2GM^rH^HC~^CH^ ■ Here, rn = 
Rs cos^(A/2) is the radius of the black hole horizon; sin A = 
q/Mh is the ratio of specific angular momentum to the 
maximal Kerr value. 

Taking SEk/SEss ~ GMnMd/A^, we obtain tn = 
{%rH /c){5Ek/SEB,e){MH /Md)\HeQ/2)\'' . Assuming a ~ 
M so ^ Rs/2, we find 



tH = 88s 



f Mh \ (Mh/MA f6Ek/SEB^ 



VIOM0 



100 



100 



(9) 

consistent with a bimodal offset of long- from short- 
duration bursts for Mj, « Mh- Note that tn/tin ~ 
Mh/Mci, as magnetic torques on the disk and black hole 
are comparable, but Jh exceeds J of the inner disk by 
a factor of the mass ratio, assuming comparable specific 
angular momenta, as for a maximal-Kerr hole. 

We remark that a black hole mass Mh of about lOM©, 
formed promptly in the hypernova scenario, is expected in 
view of the Kerr constraint Jh < GM'jj/c. For prompt 
collapse of a He-core of radius Rhci stripped of its hydro- 
gen envelope and tidally locked with a binary companion 
of mass TO (Woosley 1993; Paczyhski 1997; Iben & Tutukov 
1996) with an orbital separation ^Rhb, 



5 J \GM^3 



Here, the central density pc is that extrapolated from the 
He-mantle, with p denoting the mean density. For canon- 
ical values m = M = 2OM0, Rne = Rq and ^ = 4, 
a black hole may form promptly provided the inequal- 
ity on the left hand-side of equation (10) is satisfied; i.e. 
1 > 1.3 X 10^ {P/Pcf- This would be true, \i e.g. the He 
core is approximately fit by a polytrope with index n = 3, 
for which the Lane-Emden relationship gives pc ^ 50p 
(Kippenhahn & Weigert 1990). Thus, rapidly rotating 
black holes formed in prompt collapse tend to be massive. 

4. SUSPENDED ACCRETION AND LONG GRBS 

Long GRBs rely on the continuing presence of the torus. 
Hyperaccretion times (§2) are far shorter, however. We 
speculate that resolution of this paradox depends on the 
ability of the black hole to spin up the disk over in- 
terconnecting magnetic field-lines equivalent in poloidal 
topology to pulsar magnetospheres (Fig. 2 in van Putten 
(1999)). The black hole-torque on the torus is (adapted 
from Thorne et al. (1986)) 

1 



5t+ ^-{^h- n)sui'0{-dM^A^^ 



(11) 



for a fiux-tube of field-lines between them (counting both 
above and below the midplane). The torus now has an 
inner face interacting with the black hole, and an outer 
face interacting with infinity. The primary role of the in- 
ner interaction is to extend the lifetime of the torus by 
preventing its accretion onto the black hole. 

The total black hole torque on the torus, r+, depends 
on the associated horizon fiux 211 Ah- Approach to the 
horizon tends to increase Ah] if 211 Ah becomes large com- 
pared to the open torus flux, an equilibrium state r+ ^ r_ 
may be reached, and accretion onto the black hole stalls. 
This state is expected to be stable, since t+ decreases with 
increasing TZ: if the torus gains (loses) angular momen- 
tum, it moves outward (inward) and reduces (increases) 
T+. Such oscillations may give rise to intermittency. The 
black hole spindown proceeds on essentially the timescale 
in equation (9); the difference is that now the torus - with 
its two faces - mediates the interaction. 

More quantitatively, integration of equation (11) gives 
the net torque on the torus by the black hole as r+ « 
{Qh — ^)A'^fH/c, while the wind torque is now t_ « 
nA^f^[c\\n{eo/2)\]-^- Here, A is now (27r)-i times the 
total torus magnetic flux, and fn and /,„ are the frac- 
tions of A passing through the inner and outer light sur- 
faces, respectively, so that fn + fw ~ 1/2 — 1- The torus, 
"sandwiched" between the black hole and infinity, obeys 
an evolution equation 

2A2 r„. A 



Mdc 



1 



(12) 



(10) 



\H0o/2)\_ 

Accretion may therefore be halted - or reversed - if 
Qh :t n[l + (f^,/fHf\\n{eo/2)\-^]- Early on, fn « 1, 
and this equation will not be satisfied. Beginning near the 
midpoint of accretion when fn^fwi and provided that 
Q,H is initiated close to its maximal value, accretion may 
be suspended for an interval ^ tH^ until the value of ^h 
drops below the critical value. Since in the "suspended" 
state, the torus gains and loses almost equal quantities of 
angular momentum, significant baryonic matter may be 
carried off in the wind. 
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5. DISCUSSION AND COMPARISON WITH OBSERVATIONS 

We have revisited accretion timescales in magnetized 
black hole plus disk/torus systems, analyzing evolution 
when external torques rather than internal "turbulent" 
torques dominate. We propose two dynamical states - hy- 
pcraccretion and suspended accretion - for systems with 
low- or high-spin black holes. Around a slowly spinning 
black hole, wind torques drive the disk to a finite-time sin- 
gularity - producing short bursts. Rapidly spinning black 
holes transfer sufficient angular momentum to the disk to 
arrest matter inflow until the black hole spins down - pro- 
ducing long bursts. A bimodal duration distribution oc- 
curs, provided that « Mh- 

Our principal assumption is that external ordered 
torques dominate internal disordered ones; this requires 
significant poloidal magnetic flux. Determining how 
poloidal flux is created and when external torques prevail 
is a major unsolved problem. Numerical studies to date 
have shown the development of magnetic power spectra 
dominated by the largest scales permitted, in local disk 
simulations (e.g. Hawlcy, Gammic, & Balbus (1995)), and 
that magnetically-driven winds may develop from both 
cold and hot disks (e.g. Shibata & Uchida (1986); Stone & 
Norman (1994); Stone & Pringle (2000)), with dependence 
of the angular momentum loss rate on Eb similar to an- 
alytic predictions (Kudoh, Matsumoto, & Shibata 1998). 
Because the choice of numerical boundary conditions (par- 
ticularly for B(j,) may affect the solution at the largest scale 
(see e.g. Ustyugova ct al (1999)), very large dynamic range 
would be needed to assess the external magnetic torque. 

Our model differs from expectations for a turbulence- 
dominated disk in several ways. First, a wind-dominated 
system converts gravitational energy into Poynting flux 
rather than thermalizing it. Such a "clean" system could 



have less radiative contamination/spectral degradation 
from a low-energy thermal component. Second, a magne- 
tized wind is naturally collimated, unlike the fireball in the 
alternative scenario. Finally, only if significant poloidal 
flux is present can field geometry become unfavorable for 
matter inflow, allowing an extended disk lifetime. In tur- 
bulent disks, magnetic fields lie primarily parallel to the 
midplane, so flow onto the black hole proceeds even as 
magnetic stresses transfer spin energy and angular mo- 
mentum outward (Krolik 1999; Gammic 1999). 

Because most of the energy is liberated in the innermost 
regions, our hyperaccretion/short GRB and suspended- 
accretion/long GRB proposal implies similar intrinsic lu- 
minosities, but larger fluences in the second case. We thus 
expect short GRBs to feature afterglows similar to those 
of long GRBs, unless the environments of high- and low- 
angular momentum progenitors differ appreciably. This 
may be tested by the HETE-II mission, and is already 
supported by the afterglow detection to the 2s event GRB 
000301C (Jensen et al. 2001). 

The iron emission lines from GRB 970508 (Piro et al. 
1999) may reflect the composition of progenitor matter. 
In the suspended accretion (long GRB) state, some torus 
matter is expected to be blown to infinity by the pow- 
erful black hole torque. We associate this baryonic wind 
with line-emitting regions for long burst systems; they are 
expected to be less powerful in short GRBs. 
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